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We consider the critical behavior of the most general system of two TV-vector 

order parameters that is O(N) invariant. We show that it may a have a multi- 

£N1 i critical transition with enlarged symmetry controlled by the chiral 0(2)®0(N) 

fixed point. For N = 2, 3, 4, if the system is also invariant under the exchange 

of the two order parameters and under independent parity transformations, one 

may observe a critical transition controlled by a fixed point belonging to the 

ran model. Also in this case there is a symmetry enlargement at the transition, 

the symmetry being [SO(N) © SO(N)] ® Ci, where Ci is the symmetry group 

' of the square. 

Ph. 

Key words: N -vector model, O(N) symmetry, multicritical transitions. 



PACS: 05.70.Jk, 64.60. Fr, 75.10.Hk. 



X. 

1. Introduction 



The critical behavior of a system with a single iV-vector order parameter is well 
known [EJEJ- In this paper we investigate the critical behavior of a system with two 
N- vector parameters that is invariant under O(N) transformations and independent 
parity transformations. 

If the two iV-vector order parameters are identical, i.e. the model is symmetric 
under their exchange, the most general Landau- Ginzburg- Wilson (LGW) $ 4 Hamil- 
tonian is given by 



7i c 



+|(0 4 + ^ 4 ) + ^0V + ^(0-# 

where tp a and <fr a are two iV-dimensional vectors. This Hamiltonian is well defined 
for uq > 0, 2uq + wq > 0, and 2uq + wq + zq > 0. Hamiltonian (JTJ) is invariant under 
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the transformations 

(Z 2 )exch® (Z 2 ) par ®0(iV). (2) 

The first Z 2 group is related to the exchange transformations «-> if), while the 
second group is related to the parity transformations — > —0, ^ — > ■0, or, equiva- 
lently, — > 0, 0> — > — -0 (note that the transformation — > —0, -0 — > —-0 is already 
accounted for by the O(N) group). 

If the two order parameters are not identical and therefore the symmetry is only 

(Z 2 ) par ®0(iV), (3) 

the corresponding LGW $ 4 Hamiltonian is 

H mci = I d 3 x X - J2(9^(f> ■ + ■ djP) + ^0 2 + 

that is well defined for m > 0, vq > 0, + 2y/u v > 0, and w + z + 2^u v > 0. 
Hamiltonian (@J) has two different mass terms and thus it gives rise to a variety 
of critical and multicritical behaviors. It generalizes the multicritical Hamiltonian 
considered in Ref. [0] that has wq = and is symmetric under the larger symmetry 
group O(N) ®0(N). 

For N = 2 there is a transformation of the fields and couplings that leaves 
invariant Hamiltonians (JTJ) or (jjj). If we transform the fields as (p' a = J2b e ab4>b, 
"0a = V'a an d the couplings as 

u' = m , v'o = v , w' = w + Zq, z' = -Zq, (5) 

we reobtain Hamiltonians (0) and (j3J) expressed in terms of the primed fields and 
couplings [Hj. This implies that, for any FP with z > 0, there exist an equivalent 
one with the same stability properties and z < 0. 

Finally, if we do not require the invariance of the model under independent parity 
transformations, i.e., the model is only O(N) symmetric, we must add an additional 
quadratic term • -0 and two additional quartic terms, (0 ■ 0)0 2 and (0 ■ ifj)ip 2 . In 
this case, the general analysis becomes more complex since we have to deal with a 
multicritical theory with three quadratic terms. 

In this paper we investigate whether theories (JTJ and (j3J) have stable fixed points 
(FPs) in three dimensions. We do not determine the renormalization-group (RG) 
flow in the full theory, but rather we show that stable FPs can be identified by 
an analysis of the submodels whose RG flow is already known. We consider the 
stable FPs of the submodels and determine their stability properties with respect 
to the perturbations that are present in the complete theory. In this way we are 
able to identify three stable FPs. For any N, there is an 0(2) ® O(N) symmetric 
FP. This FP may be the relevant one for models with z > 0, which may therefore 
show a symmetry enlargement at the (multi) critical transition. For N = 2 there is 
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an equivalent 0(2) £g> 0(2) symmetric FP with z < 0, a consequence of symmetry 
(jSJ), which may be the relevant one for models with zq < 0. For iV = 2,3,4 we 
find that 7i cr — but not the multicritical theory 7i mcr — has another stable FP that 
belongs to the so-called mn model [ |3] with n = 2 and m = N. Also in this case 
there is a symmetry enlargement at the transition: the FP is symmetric under the 
group [SO(N) © SO(N)] ® C 2 where C 2 is the symmetry group of the square. It is 
interesting to note that the chiral 0(2) ®0(N) FP is also stable if we do not require 
the model to be invariant under independent parity transformations. Indeed, the 
additional terms (0 ■ ip)<j) 2 and (0 • ip)ip 2 are irrelevant perturbations at the chiral 
FP. 

In the analysis we mainly use the minimal-subtraction (MS) scheme without 
e expansion (henceforth indicated as 3d-MS scheme) in which no e expansion is 
performed and e is set to the physical value e = 1 [ 6 . In order to generate the 
relevant perturbative series we use a symbolic manipulation program that generates 
the diagrams and computes symmetry and group factors. For the Feynman integrals 
we use the results reported in Ref . [ Ej . In this way we obtained five-loop 3<i-MS 
expansions. 

2. Mean-field analysis 

The mean-field analysis of the critical behavior of Hamiltonian 7i cr is quite 
straightforward. If r > the system is paramagnetic, with = -0 = 0. For r < 
there are three possible low-temperature phases: 

(a) For wo > 2uo and z > 2uo — wo, we have ^ and ip — (or viceversa). 
The corresponding symmetry-breaking pattern is (Z 2 ) CX ch <8> (Z 2 ) par <8> O(N) — > 
(Z 2 ) pax ®0{N-l). 

(b) For z < and — 2m < wo + z < 2w , we have </> = ^ 0. The corresponding 
symmetry-breaking pattern is (Z 2 ) exch Cg>(Z 2 ) pai .®0(iV) — > (Z 2 ) exch ®0(iV — 1). 

(c) For z > and — 2w < wq < 2m , we have \<p\ = \ip\ ^ 0, • ip = 0. 
The corresponding symmetry-breaking pattern is (Z 2 ) cxc h ® (Z 2 ) par ® O(N) — > 
(Z 2 ) cxch ®0(iV-l). 

The analysis of the mean-field behavior of 7Y mcr is presented for N = 4 in Ref. [|H] 
and it is easily extended to the present case. There are three possible phase diagrams: 

(al) For zq < and — 2^/uqVq < w + z < 2 v /w t> , the multicritical point is 
tetracritical, see Fig. ^ Phase 1 is paramagnetic with = ip = 0, in phase 2 
7^ and ip = 0, while in phase 3 the opposite holds, = and ip ^ 0; in 
phase 4 ^ 0, -0 with (f)\\ip. All transitions are of second order. Transitions 
1-2 and 1-3 are associated with the symmetry breaking Z 2 <g> O(N) — » Z 2 <g> 
0(N — 1). In the presence of fluctuations these transitions belong to the O(N) 
universality class. Transitions 2-4 and 3-4 are associated with the symmetry- 
breaking pattern Z 2 eg) 0(N — 1) — > 0(N — 1). In the presence of fluctuations 
they should belong to the Ising universality class. 
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Figure 1. Possible multicritical phase diagrams. Thin lines indicate second-order 
transitions, while the thick line in case (b) corresponds to a first-order transition. 
"Is" indicates an Ising transition. 



(a2) For z > and — 2-^/wnfn < u>o < Zy/uoVo the multicritical point is tetracritical, 
see Fig. [H Phases 1,2, and 3 as well as transitions 1-2 and 1-3 are identical 
to those discussed in case (al). In phase 4 ^ 0, if) ^ with • if) = 0. All 
transitions are second-order ones. Transitions 2-4 and 3-4 are associated with 
the symmetry-breaking pattern Z 2 (g> 0(iV — 1) — > Z 2 <8> 0(iV — 2) and, in the 
presence of fluctuations, they should belong to the 0(N — 1) universality class. 

(b) For w Q > 2y/u Q vo and wo + zo > 2y/u v the multicritical point is bicritical, 
see Fig. H Phases 1, 2, and 3 as well as transitions 1-2 and 1-3 are identical to 
those discussed in case (al). The transition between phases 2 and 3 is of first 
order. 



3. Analysis of some particular cases 

3.1. Particular models and fixed points 

The three-dimensional properties of the RG flow are determined by its FPs. 
Some of them can be identified by considering particular cases in which some of the 
quartic parameters vanish. For H cr we can easily recognize two submodels: 
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(a) The 0(2) ® O(N) model with Hamiltonian [Ej 



n 



ch 



d d x 



4! V Z_^ « ; ^ 4! 



l,j a 



(6) 



where $ a j is an A/" x 2 matrix, i.e., a = 1, . . . , N and i — 1,2. Hamiltonian (|T|) 
reduces to © for 2ii -wo--2o = 0, if we set $ al = <f) a , $ a2 = ip a , u = #i,o, w = 
2(g li0 —g 2) o), and z = 2g 20 . The properties of 0(2)®0(N) models are reviewed 
in Refs. [EH 121 El EI] - I n three dimensions perturbative calculations within the 
three-dimensional massive zero- momentum (MZM) scheme [E1E] an< ^ within 
the 3d-MS scheme [ E] indicate the presence of a stable FP with attraction 
domain in the region g 2 ,o > for all values of N (only for N = 6 the evidence is 
less clear: a FP is identified in the 3<i-MS scheme but not in the MZM scheme). 
For N = 2, these conclusions have been recently confirmed by a Monte Carlo 
calculation [E|- On the other hand, near four dimensions, a stable FP is found 
only for large values of N, i.e., N > N c = 21.80 - 23.43e + 7.09e 2 + 0(e 3 ) [ 
13 El El E] • A stable FP with attraction domain in the region g 2 ,o < exists 
for N = 2 (it belongs to the XY universality class) [E|, for N = 3 (Ref. [CH]), 
and N = 4 (Ref. [HJ). Note that nonperturbative approximate RG calculations 
have so far found no evidence of stable FPs for N = 2 and 3 [El El 120] ■ in 
the following we call the FP with g 2 > chiral FP (we indicate it with gl ch , 
92 ch > 0)> while the FP with g 2 < is named collinear FP and indicated with 

(b) The so-called mn model with Hamiltonian [|3] 




4! 1 



(7) 



where $ a j is an m x n matrix, i.e., a 



m and z 



AT, and z = 0, if we set $ 



al 



n. Hamiltonian 

= 0a, $a2 = Ipa, 



(HJ) reduces to ([0]) for n = 2, m 

u o — gi,o + #2,0 an d w = 2gi Q. A stable FP is the 0{m) FP with g x = and 
#2 = 9mi where ^ is the FP value of the renormalized coupling in the 0(m) 
model. In App. El we show that the model has a second stable FP with g 2 < 
for n = 2 and m — 2, 3, and 4. We name this FP the mn FP and we label the 
corresponding coordinates by g* mn and g 2mn . 

The presence of these two submodels that have one parameter less than the original 
model imply that the quartic parameter space splits into four regions such that the 
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RG flow does not cross the two planes z = and 2u — w — z = 0. Note that, for 
N = 2, because of symmetry we should only consider the region z > 0. 

The results for models (a) and (b) allow us to identify four possible FPs that are 
candidates for being stable FPs of the full theory: 

(1) u = g* hch , w = 2(g* lch - #* ch ), z = 2g* ch ; this FP may be the stable FP of the 
trajectories that start in the region zq > 0. 

(2) u = g* cl , w = 2(g* d — <?2ci)> z = cii this FP may be the stable FP of the 
trajectories that start in the region z < 0. 

(3) u = g* N , w = 0, z = 0; this FP may be the stable FP of the trajectories that 
start in the region 2uq — wq — z$ > 0. 

(4) u = g\ mn + g* 2) mn , w = 2g\ mn , z = for iV = 2, 3, 4; this FP may be the stable 
FP of the trajectories that start in the region 2uq — wq — zq < 0. 

Note that, because of symmetry for N = 2 there is also a chiral (resp. collinear) 
FP with z < (resp. z > 0). 

The analysis of the particular cases of Hamiltonian (|IJ) is very similar. There are 
two relevant submodels: 

(a) The forementioned 0(2) ® O(N) model for u = v and w + z = 2u . The 
identication is u = v = g lfi , w = 2g lfi - 2g 2>0 , and z = 2g 2fi . 

(b) The O(N) © 0(N) model [0: 

= I d 3 X ^ ^(^0 • + ■ drf) + j^ 2 + + 

. fl,0 ,4 , ^2,0 ,4 , ^3,0 ,2 i 2 /q\ 

+ir + ir^ + ^r ^ J • (8) 

Hamiltonian 7i mcr reduces to this model for zq = 0, with the obvious identifica- 
tion of the parameters. Hamiltonian (JHJ) describes the multicritical behavior of 
a model with two iV-vector order parameters that is symmetric under indepen- 
dent O(N) transformations of the two order parameters, i.e. that is invariant 
under the symmetry group O(N) © O(N) [E]. In the case we are interested 
in, i.e. for N > 2, the stable FP is the decoupled FP [HE1I221, i.e., h = 0, 
h = h = 9*n ( see also A PP- EJ- 

Note that in this case the 0(2) © O(N) model has two parameters less than the 
original one and thus its presence does not imply any separation of the RG flow. 
Instead, the second model implies that the quartic parameter space splits into two 
regions such that the RG flow does not cross the plane z = 0. The analysis of the 
possible FPs is identical to that presented above, since the FPs we have identified 
are exactly those we have already described. 
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3.2. Stability of the 0(2)®0(iV) fixed points 

In this section we study the stability properties of the two FPs that appear in 
the 0(2) ® O(N) model. For this purpose we need to classify the perturbations of 
the 0(2) (g> 0(N) model that do not break the 0(N) invariance. The multicritical 
Hamiltonian (J3J) can be rewritten as 

n mcr = n ch + \r 2 , 2 v^ + \h A v^ + lhiv^\ (9) 

where 7i c h is the 0(2) ® 0(iV)-symmetric Hamiltonian (JBJ) with r = (r x + r 2 )/2, 
#1,0 = (2u + 2v + w + z )/6, 02,0 = (2u + 2v -2w + z )/6, and r 2 , 2 = (ri-r 2 )/2, 
/ 4>2 = ( Uo - Wo )/2, and / 4 , 4 = (u + v -w -z )/6. Here, 1/( 2 ' 2 \ 4 ' 4 ), and V^< 2 \ are 
respectively a quadratic term that transforms as a spin-2 operator under the 0(2) 
group and two quartic terms that transform as a spin-4 and as a spin-2 operator 
respectively. Their explicit expressions are: 

V^ 4 ' 4 ) = (0 2 ) 2 + (t/; 2 ) 2 - 20 2 ^ 2 - 4(0 • VO 2 . (10) 

A detailed description of all possible perturbations of the 0(2) ®0(N) FP the leave 
invariant the 0(N) group can be found in App. B of Ref. [0]. Note that for 7i cv we 
have r 2j 2 = fi,2 — 0, so that one must only consider the spin-4 quartic perturbation. 

Let us first discuss the chiral FP (a) that has g\ ch > 0. In order to estimate the 
RG dimensions 4j2 and of the above-reported perturbations, we computed the 
corresponding five-loop MS series and we analyzed them within the 3c?-MS scheme. 

The perturbative series, that are not reported here but are available on request, 
were analyzed using the conformal-mapping method and the Pade-Borel method, 
following closely Ref. [123; to which we refer for details. The error on the conformal- 
method results takes into account the spread of the results as the parameters a and 
b are varied (cf. Ref. for definitions) and the error due to the uncertainty of the 
FP location (we use the estimates reported in Refs. [H2JE31II21I11)- The results of the 
analyses using the conformal-mapping method are reported in Table ^ Completely 
consistent results are obtained by using Pade-Borel approximants. As it can be seen, 
04 >2 and 7/4,4 are always negative, indicating that the chiral FP is stable for any N. 
This FP is therefore expected to be the relevant FP whenever the RG flow starts in 
the region z > 0. For N = 2, symmetry (JHJ) implies that a chiral FP [the equivalent 
one that is obtained by using ©] may also be reached from the region z < 0. 

It is of interest to compute also the RG dimension t/ 2i2 of the quadratic pertur- 
bation. For the multicritical Hamiltonian it is related to the crossover exponent <p\ 
cj) = z^/2,2) where v is the correlation- length exponent at the chiral FP (see Refs. [ 
O CH U2] for numerical estimates). The exponent t/ 2j2 has already been computed 
for several values of iV in Ref. indeed, 2/2,2 = 04; where 04 is the RG dimension 
of the operator O^ defined in App. C of Ref. [Hj. Numerical estimates, taken from 
Ref. [Hj, are reported in Tabled] 
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2/2,2 


2/4,2 


2/4,4 


ch,2 


1.34(15) 


-1.6(1.1) 


-0.9(4) 


ch,3 


1.21(9) 


-1.4(8) 


-1.0(3) 


ch,4 


1.17(8) 


-1.1(5) 


-1.0(3) 


ch,6 


1.13(9) 


-0.9(4) 


-0.9(2) 


ch,8 


1.13(8) 


-0.9(5) 


-0.9(2) 


ch,16 


1.08(2) 


-1.0(3) 


-0.94(11) 


ch,oo 


1 


-1 


-1 


cl,2 


1.9620(8) 





0.532(12) 


cl,3 


2.05(15) 




0.9(3) 


cl,4 


2.05(15) 




0.9(7) 



Table 1. Estimates of the RG dimensions 7/2,2, 2/4,2j and 2/4,4 of the operators 
V& 2 \ V^ 2 \ and at the chiral (ch) FP and at the collinear (cl) FP. They 

have been obtained from a conformal-mapping analysis of the corresponding 3d- 
MS 5-loop perturbative expansions. The results for 7/2,2 are taken from Ref. [ 
0]- The results at the collinear FP for N = 2 have been computed by using the 
mapping with the XY model and the results of Ref. [ |^| . 



Now, let us consider the collinear FP (b) that has g% C \{N) < for 2 < N < 4. For 
N = 2 the RG dimensions at the collinear FP can be related to the RG dimensions 
of operators in the XY model. Indeed, the 0(2)® 0(2) collinear FP is equivalent to 
an XY FP. The mapping is the following. One defines two fields a, and bi, i = 1,2, 
and considers [E] 

0n = (ai - b 2 )/V2, 
022 = (ai + b 2 )/V2, 
012 = (h - a 2 )/V2, 

02i = (61 + a 2 )/V2. (11) 

At the collinear FP, fields a and b represent two independent XY fields. Using 
this mapping it is easy to show that V^ 4,2 ^ ~ Of ,1 \a)bj, Of (b)aj, and V^ 4,4 ^ ~ 
T n (a)T u (b), T 12 (a)T 12 (b), where 

®T ' l \ a ) = a i a2 i T i:j (a) = aidj - ^Sija 2 . (12) 

(3 1) 

Thus, if ys t i and y 2 are the RG dimensions of 0\ ' and Ty in the XY model, we 
have 

2/4,2 = Vh + 2/3,1 - 3, 2/4,4 = 2?/2 - 3. (13) 

By using the equations of motion, one can relate 0\ to [ Elj and obtain y 31 = 
3 — yh, so that 7/4 2 = exactly (this holds in three dimensions; in generic dimension 
d, 2/4,2 = 3 — d). For y 2 we can use the result reported in Ref. [ |2H, obtaining 
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= 0.532(12). The analysis of the perturbative series gives results that are fully 
consistent: 2/4,2 = 0.0(1), 2/4,4 = 0.57(4). 

In order to determine 2/4,2 and 2/4,4 for N = 3 and 4 we analyzed the corresponding 
5-loop 3rf-MS expansions. The results for 2/4,4 are reported in Tabled They indicate 
that 2/4,4 is positive, which implies that the spin-4 quartic perturbation is relevant 
and therefore the collinear FP is unstable. We do not quote any result for 2/ 4 2 . The 
perturbative analysis does not allow us to obtain any reliable result: the estimates 
vary significantly with the parameters b and a and with the perturbative order. 

It is interesting to observe that, on the basis of the group-theoretical analysis 
reported in App. B of Ref. [Hj, any quartic perturbation of the chiral FP that leaves 
invariant the O(N) symmetry is a combination of spin-2 and spin-4 operators. Thus, 
the results presented here indicate that the chiral 0(2) ® O(N) FP is stable under 
any perturbation that preserves the O(N) symmetry. In particular, it is also stable 
under a perturbation of the form (<j)-if))(a<f) 2 + bip 2 ) that may arise if the model is not 
invariant under independent parity transformations. Indeed, such a term is nothing 
but a particular combination of spin-2 and spin-4 perturbations. In the notations of 
App. B of Ref. [E] (note that M and N of Ref. [0] should be replaced by iV and 2 
respectively) we have 

(0 • 2/;)(a0 2 + 62A 2 ) = i(a + b)0^ l) + i(a - b)0$$. (14) 

Moreover, the additional quadratic term <f> ■ ip is nothing but a component of the 
spin-2 quadratic term that breaks the 0(2) group, so that the associated crossover 
exponent is again = z/2/2.2, 1/2,2 being reported in Tabled This is a general re- 
sult that follows from the analysis of Ref. [0j: any quadratic perturbation of the 
0(2) (gi O(N) FP that does not break the O(N) invariance is a combination of 
the components of the spin-2 quadratic operator. Thus, any perturbation is always 
associated with the same crossover exponent (f> = ^2/2,2- 

3.3. Stability of the decoupled O(N) © O(N) fixed point 

We now consider FP (3) discussed in Sec. 13.11 In order to check its stability, we 
must determine the RG dimensions at the FP of the perturbations 

P E = </>V, Pt = J2 T *M T 1>w ( 15 ) 

u 

where T^j = 4>i<Pj — ^Sij4> 2 . Simple RG arguments show that the RG dimensions 
are given by 

2 On 

Ve = 3 = — , y T = 22/2 - 3, (16) 

u N v N 

where and are the critical exponents of the 3-dimensional O(N) universality 
class (see Ref. [|2] for a comprehensive review of results), while 2/2 is the exponent 
associated with the quadratic spin-2 perturbation in the 0(N) model [ |2Hl EH I2H] ■ 
Since < for iV > 2 we have 2/^ < 0, i.e. the perturbation P E is always 
irrelevant. As for yx we can use the results reported in Refs. [ 123 I2HI- The 
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spin-2 exponent is equal to y 2 = 1.766(6), 1.790(3), 1.813(6) for AT = 2, 3, 4 and 
increases towards 2 as N — > oo. Correspondingly jjt = 0.532(12), 0.580(6), 0.626(12), 
increasing towards 1 as iV — > oo. It follows that Pt is always relevant. Thus, the 
decoupled FP is always irrelevant. 

3.4. Stability of the mn fixed point 

Here, we wish to consider the stability of the mn FP. For this purpose we must 
consider the two perturbations 

Pi = (0 2 ) 2 - (^ 2 ) 2 , P2 = (0 • ^) 2 - ^V- (17) 

Note that Pi is not symmetric under interchange of <fi and ip and is therefore not 
of interest for 7i cr . The corresponding RG dimensions y\ and y 2 are computed in 
App-El yi = 0.4(3), 0.2(2), 0.2(2) for AT = 2, 3, 4; y 2 = -0.9(5), -1.0(8), -0.8(5) 
for the same values of N. They indicate that P\ is relevant and P 2 is irrelevant at the 
mn FP. Therefore, the mn FP is a stable FP for T~C cr (only P 2 should be considered 
in this case) and an unstable one for Ti. mcT - 

4. Conclusions 

In this paper we have investigated the critical behavior of systems described by 
Hamiltonians and (pEfl. We find that 7i cr has three possible stable FPs: for any N, 
except possibly N = 6 (for such a value of A" the evidence of this FP is less robust 
[112]). there is the 0(2) ® O(N) chiral FP that is relevant for systems with zq > 0; 
for N = 2 there is a stable chiral FP with z < [equivalent to the previous one by 
symmetry (JSJ)], that is relevant for systems with z < 0; for iV = 2, 3, 4, there is the 
mn FP that is relevant for systems with 2uq — Wq — z < 0. In the multicritical theory 
P)) only the chiral FPs are stable. Thus, systems with z > (or, for N = 2, with 
Zq 7^ 0) may show a multicritical continuous transition with the larger 0(2) ®0(N) 
symmetry. 

It is interesting to note that the most general 0(iV)-invariant LGW Hamiltonian 
for two iV-vector parameters includes other couplings beside those present in (0) 
and (jlj. One should consider 

Hcr.ext = K r + y0^ + ^(0^)(0 2 + ^ 2 ), (18) 
^mcr.ext = ^mcr + y <P ' 4> + (4> ' ^) ( a 10 2 + «2^ 2 ) , (19) 

depending whether one wants to preserve the symmetry under the exchange of the 
two fields. As we discussed in Sec. 13.2] the chiral FP is a stable FP also for these 
two extended models. 

Hamiltonians and (TTS|) have two mass parameters and thus symmetry en- 
largement can be observed only at the multicritical point, where the singular part 
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of the free energy has the form 

i^sing = iff (20) 

where /i 4 and fi g are two linear scaling fields (linear combinations of the temperature 
and of another relevant parameter), a and <fi are the specific- heat and the crossover 
exponents at the 0(2) ®0(N) model. Note that the same expression, with the same 
a and <f), applies to both models, apart from nonuniversal normalization constants. 
In 7i mcr ,ext there are three quadratic parameters and thus the chiral multicritical 
point can be observed only if three relevant parameters are properly tuned. The 
singular part of the free energy becomes 

-Fsing = /i(" a /(/i g l/i 4 , fJ-gZfJ-t )> ( 21 ) 

where [L g \ and \i g2 are two linear scaling fields associated with the same crossover 
exponent 0. 

A. The mn model: new fixed points 

In this Appendix we consider the mn model defined by Hamiltonian (J7j). focusing 
on the case n = 2 that is of interest for the present paper. Within the e expansion one 
finds four FPs, the stable one being the 0{m) FP with g x = and g 2 = g^, where 
g* m is the FP value of the renormalized zero- momentum coupling in the 0(m) model, 
see Refs. [ [27] and references therein. For m = 2 (and n = 2) the mn model is 
equivalent to the 0(2) ®0(2) model defined by Hamiltonian ©. For this model, 
the results of Refs. [H31II2] indicate the presence of a new FP that is not predicted 
by the e-expansion analysis. Because of the mapping, this implies the presence of 
a new FP in the mn model with g 2 < 0. In the MZM scheme the results of Ref. [ 
13j imply the presence of a stable FP at g\ = 4.4(2) and g 2 = —4.5(2), where the 
renormalized couplings g± and g 2 are normalized so that g± = 3gifi/(lQ7rR 2m m), g 2 = 
3^2,0/ {lQ7[R m m) at tree level (m is the renormalized zero-momentum mass), where 
Rk = 9/(8 + k). In the 3rf-MS scheme, by using the results of Ref. [H2], we obtain 
gx = 2.25(13) and g 2 = -2.31(21), where g { = g^-' 7 1 A d with A d = 2 d " V/ 2 r(d/2). 
It is thus of interest to check whether additional FPs are also present for other values 
of m > 2. As we shall show below we find an additional FP for m = 3 and m = 4. 
For m > 5 no new FP is found. 

In order to check for the presence of additional FPs we considered the six-loop 
MZM expansions of Ref. [ |2H] and we generated 5-loop 3<i-MS expansions. For the 
analysis we used the conformal-mapping method: the position of the Borel singu- 
larity in the MZM scheme is reported in Ref. [I^Eh while in the 3c?-MS we used its 
trivial generalization. The two (3 functions were resummed by using several different 
approximants depending on two parameters, b and a (see Ref. [123] for definitions). 
For simplicity, each time we resummed the two (3 functions by using the same b 
and a and then determined their common zeroes. In principle, it would have been 
more natural to consider different values of b and a for the two (3 functions and all 
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Table 2. Results for the mn model for n = 2 in two different schemes. The index 
in column "scheme", 41, 51, 61, refers to the number of loops. We report the 
coordinate of the FP g\, g2, the percentage of approximants that find the zero 
(pfp)j and the percentage of approximants that indicate that the FP is stable 
(p s t). In the column "info" we report the number of approximants that give real 
(first number) and complex eigenvalues (second number). 

possible combinations. However, as we already tested in the 0(2) <S> O(N) model, 
the two choices give fully equivalent results. In the analysis we used a = —1, 0, 1, 2 
and b = 4, 6, . . . , 14, which appeared to be an optimal choice. We report the results 
in Table El For comparison, we also performed the analysis for m = 2, obtaining 
results completely consistent with those reported above. In the table we also give the 
percentage of cases in which a FP was found (pfp) and in which this FP was stable 
(p st ). Finally, we also indicate the number of cases in which the stability eigenvalues 
were real or complex. 

For m = 2 and m = 3 the presence of a new FP is unambiguous. Essentially 
all considered approximants at five and six loops in both schemes give a stable FP. 
For m = 4, the percentages are smaller, although the overall results are still in 
favor of a new stable FP. For m > 5 there is essentially no evidence. As far as the 
stability eigenvalues, for m = 2 they are complex, in agreement with the results for 
the 0(2) ® 0(2) model [E3 EH 02] ■ F° r m — 3 and m = 4 the numerical results 
favor real eigenvalues instead. It is interesting to note that the 3c/-MS FPs lie at the 
boundary of the region in which the expansions are Borel summable, g\+g2 > 0. This 
is not the case for the MZM ones (Borel summability requires R2m9i + Rm92 > 0). 
Thus, the 3<i-MS results should be more reliable in these models. 

The mn model for n = 2 is invariant under the group [SO(m) © SO(m)] <S> 
C2 where C2 is the symmetry group of the square. We now consider two quartic 
operators that break such a symmetry: 

P 1 EE ($! • $!) 2 - ($ 2 • $ 2 ) 2 , (22) 
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P 2 = ($ 1 -$ 2 ) 2 --($ 1 -$ 1 )($ 2 -$ 2 ), (23) 
m 

where the scalar products are taken in the 0(m) space. The first operator is the 
only quartic one that preserves the continuous symmetry and breaks C 2 — > Z 2 ©Z 2 , 
while the second preserves C 2 but breaks SO{m) © SO{m) — > SO{m). Note that 
in general Pi mixes with the lower-dimensional operator $! • $ x — $ 2 • $ 2 . Such a 
mixing should be taken into account in the MZM scheme, but does not occur in the 
massless MS scheme. The operators Pi and P 2 are relevant for the stability of the 
FPs of the mn theory in larger models with smaller symmetry group. 

We computed the anomalous dimensions of Pi and P 2 at the new FPs of Tabled 
by analyzing the corresponding 5-loop 3d-MS series. The exponent y\ was obtained 
from the analysis of the inverse series 1/yi; the direct analysis of the series of y\ was 
very unstable. For j/ 2 we used instead the corresponding series. The results were not 
very stable and should be taken with caution. They are: 

m = 2: yi = 0.4(3), j/ 2 = -0.9(5); 

m = 3: y x = 0.2(2), y 2 = -1.0(8); 

m = 4: yi = 0.2(2), ?/ 2 = -0.8(5). 
It is interesting to note that the mn model is a submodel of the multicritical 
Hamiltonian © for n = r 2 and / 1>0 = / 3 ,q if we set $ al = a , $ a2 = / 1(Q = 
/3,o = ^1,0 + ^2,0; an d / 2 ,o = 2gi,o- This implies that the new FPs may be relevant for 
the multicritical behavior of 7i mcr)2 . To investigate this possibility we must compute 
the anomalous dimension of the operator that breaks [SO(m) © SO{m)} © C 2 — > 
0{m) ffiO(m), i.e., the operator P x . As it can be seen, y x > in all cases, indicating 
that the mn FP is unstable in the full theory. Thus, the decoupled FP appears to 
be the only stable FP of the multicritical model @ [ 1221 12H • 
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